Vorticity Knot in Two-component Bose-Einstein Condensates 
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We demonstrate the existence of the hehcal vortex solution in two-component Bose-Einstein 
condensates which can be identified as a twisted vorticity flux. Based on this we argue that the 
recently proposed knot in two-component Bose-Einstein condensates can be interpreted as a vorticity 
knot, a vortex ring made of the helical vortex. This picture shows that the knot is made of two 
quantized vorticity fluxes linked together, whose topology 7r3(S^) is fixed by the linking number 
of two vorticity fluxes. Due to the helical structure the knot has both topological and dynamical 
stability. We estimate the energy of the lightest knot to be about 3 x 10^'^ eV. 
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The topological objects, in particular finite energy 
topological objects, have played an important role in 
physics 0, 01 • In Bose-Einstein condensates (BEC) the 
best known topological objects are the vortices, which 
have been widely studied in the literature. Theoretically 
these vortices have successfully been described by the 
Gross-Pitaevskii Lagrangian. On the other hand, the re- 
cent advent of multi-component BEC (in particular the 
spin- 1/2 condensate of *^Rb atoms) has widely opened a 
new opportunity for us to study novel topological objects 
which can not be realized in ordinary (one-component) 
BEC H, 01 ■ This is because the multi-component EEC 
naturally allows a non-Abelian structure which acco- 
modates a non-trivial topological objects, in particular 
a topolgical knot which is very similar to the knot in 
Skyrme theory |^ |^ . 

Indeed recently many authors have proposed the ex- 
istence of a knot in Gross-Pitaevskii theory of two- 
component BEC |3)Hls!l- The purpose of this report is to 
show that this knot is nothing but a vorticity knot which 
is made of two vorticity fluxes linked together. Further- 
more, we show that the knot is topological, whose topology 
713(6'^) is fixed by the Chern-Simon index of the velocity 
potential of the condensate. To show this we first present 
a helical vortex solution in two-component BEC which is 
periodic in z-coordinate, and construct a helical vortex 
ring by bending it and smoothly connecting two periodic 
ends together. We show that this vortex ring becomes 
the vorticity knot whose quantum number is fixed by the 
Chern-Simon index of the velocity potential, which de- 
scribes the linking number of two vorticity fluxes. 

This picture tells that the knot has both topological 
and dynamical stability. The topological stability follows 
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from the fact that two linked vorticity fluxes can not be 
disconnected by any smooth deformation of the field con- 
figuration. The dynamical stability follows from the fact 
that the knot necessarily has a net velocity flux along 
the knot, and thus a non-vanishing angular momentum 
around the knot. This creates a repulsive stablizing force 
against the collapse of the knot. This provides the dy- 
namical stability of the knot. 

The knot that we discuss here are very similar to the 
knot in Skyrme theory 0, 0, ^ . Just as the knot in 
Skyrme theory is a vortex ring made of the helical mag- 
netic vortex, our knot here is a vortex ring made of the 
helical vorticity vortex. So it is crucial that we have the 
helical vortex to demonstrate the existence of the vortic- 
ity knot in two-component BEC. 

To construct the desired vortex solution let the two- 
component BEC be a complex doublet (j) = ((j}i, 02 ), and 
consider the Lagrangian 

C = i-cfi^dtcj) - -^\d^'t>? + Mi'/'I^i + M24</'2 

-^('^I'/'i)' - Ai2(40i)(4'^2) - ^(402)^ (1) 

where /i^ are the quadratic coupling constants and \ij 
are the quartic coupling constants which are determined 
by the scattering lengths Oij 

This is an obvious generalization of one-component 
Gross-Pitaevskii Lagrangian to the two-component BEC. 
Notice that here we have neglected the trapping poten- 
tial, because we are assuming that the range of the trap- 
ping potential is much larger than the size of tpological 
objects we are interested in. 
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Clearly the Lagrangian has a global U{1) x U{1) sym- 
metry. But one could simplify it because experimen- 
tally the scattering lengths often have almost the same 
value. For example, for the spin 1/2 condensate of ®^Rb 
atoms, all aii are about 5.5 nm and differ by only about 
3 % or so 0, 0|. In this case one may safely assume 
All — Ai2 — A22 — A. With this the Lagrangian is writ- 
ten as 



-5^402, 



Mn2 

A'' 



(3) 



where fj, = fii and Sfj, = fJ-i — fJ-2- Notice that the La- 
grangian has a global U (2) symmetry when Sfj, = 0. So 
the Sfj, interaction is the symmetry breaking term which 
breaks the global U{2) symmetry to U{1) x U{1). This 
means that even when Sfi ^ the Lagrangian has an ap- 
proximate U{2) symmetry. Physically 6fi can be viewed 
to represent the difference of the chemical potentials be- 
tween (jji and (f>2, so that it does not vanish when the 
chemical potentials are different. 



With 



V2 



pC, iCH = 1) 



(4) 



the Lagrangian Q gives the following Hamiltonian in the 
static limit (in the natural unit c = h = 1), 
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where A = 4]VP\, S^i'^ = 2M6fj,, pi = 4fiM/\, and we 
have normalized p to {^/2M /K)p. The Hamiltonian ((SJ 
can be expressed as 



n 



Sp 



C2C2}, 



(6) 



where p = p/po and di = di/VXpo. This tells that the 
physical unit of the Hamiltonian is Apg, and the physical 
scale K of the coordinates is l/VApo- This is comparable 
to the correlation length ^ — l/y/2pM. Indeed we have 
K = 1/^/2. 

From the Hamiltonian we have 
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FIG. 1: The helical vortex in the Gross-Pitaevskii theory of 
two-component BEG. Here we have put m = 1, m' = — 1, n = 
l,n' = 0, fc = 0.25/ac, and g is in the unit of k. Dashed and 
solid lines correspond to 5p/p = 0, and 0.1 respectively. 



To obtain the vortex solution, we choose the ansatz 

p = p{q)-, 

C = exp(— ^7) ^, 7 = n'lf + m'kz, 
\ sm ■ 



■ exp(— zn(p — imkz) ' 
fiQ) 



(8) 



Now, with n' — Q and m' 



P+-P 

Q 



is reduced to 



7/^ + ^ 



_ m^e Sp') sin^ l)p= _ pDp^ 



A 

2' 



Q Pj \Q 
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,2 " 
= 0. 



rri^k^ — Sp'' I sin / 



(9) 



So with the boundary condition 

p'(0) - 0, p{^) = po, /(O) = TT, /(^) = 0, (10) 

we can solve @. With m = n = 1 we obtain the twisted 
vortex solution shown in Fig. ^ 

The untwisted non-Abelian vortex solution has been 
discussed before ,9j , but the twisted vortex solution here 
is new. Notice that when Sp^ — 0, there is no untwisted 
vortex solution because in this case the vortex size be- 
come infinite. But remarkably the helical vortex exists 
even when Sp^ = 0. This is because the twisitng reduces 
the size of vortex tube. 

In Skyrme theory the helical vortex is interpreted as 
a twisted magnetic vortex whose flux is quantized 0, ^| 
Now we show that the above vortex is a twisted vorticity 
vortex. To see this notice that the non-Abelian structure 
of the vortex is represented by the doublet Moreover, 
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FIG. 2: The supercurrent (in one period section in z- 
coordinate) and corresponding magnetic field circulat- 
ing around the cylinder of radius q of the helical vortex in 
two-component BEG. Here m = l,m' = — l,n = l,n' = 0, 
k = 0.25/k, and q is in the unit of k. The current density jip 
is represented by the dotted line. 



FIG. 3: The supercurrent if and corresponding magnetic field 
flowing through the disk of radius g of the helical vortex 
in two-component BEG. Here m = l,m' = — l,n = l,n' — 0, 
k = 0.25/k, and g is in the unit of k. The current density jz 
is represented by the dotted line. 



the velocity field of the doublet is given by 

= ^(cos/(£i) + l){ndi,'f + mkd^z) + 9^7, (11) 
which generates the vorticity 

+mk{d^QdyZ - di,gdfj_z) 
This has two vorticity fluxes, (j)z along the z-axis 
4'z ^ / Vg^gdgdip = 27rn, 



(12) 



(13) 



and around the the z-axis (in one period section from 
^ = to z = 27r//c) 



27r/fc 



Vzgdgdz = — 27rm. 



(14) 



Obviously they are quantized. As importantly they are 
linked together, and have the linking number mn. 

Furthermore, just as in Skyrme theory, these fluxes 
can be viewed to originate from the helical supercurrent 
which confines them with a built-in Meissner effect 



sin/ 



COS/y2 



sin / ' 
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+mk{f+--^f + -f)d, 
smj g 



(15) 



This produces the supercurrents (in one period section 
from z = to z = 2n/k) around the z-axis 



27rn sin / • 

i0 = / 



^— oo 



and iz along the z-axis 



iz — 2Trmkgf sin / 



(16) 



(17) 



The vorticity fluxes and the corresponding supercurrents 
are shown in Fig. [3 and Fig. El This is strikingly sim- 
ilar to what we find in the magnetic vortex in Skyrme 
theory 11]. This tells that the helical vortex is nothing 
but the twisted vorticity fiux confined along the z-axis 
by the velocity current, whose flux is quantized due to 
the topological reason. We emphasize that this interpre- 
tation holds even when the is not zero. 

We can estimate the energy of the helical vortex. For 
®^Rb we have 

M ~ 8.1 X IQi" ey, A ~ 1.68 X 10"^ (nTO)2, 

^ ~ 3.3 X 10~i2 eV, Sfj. ~ 0.1 /i. (18) 

So, with m = n = 1, to' = —1, rt' = and k — 0.25/k, 
we find numerically that the energy per one periodic sec- 
tion (from z = to z = 2^ /k) is given by 



E ~ 270.987 



2.29 X 10-3 eV. 



1492.4 po 



(19) 



As we will see later, the lightest knot could have an en- 
ergy comparable to this energy. 

Notice that the vorticity (|12|l is completely fixed by 
the CP^ field ^, because it does not depend on the U{1) 
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phase 7 of C- Moreover ^ naturally defines a mapping 
from the compactified xy-plane to the target space 
S^. This means that our vortex has exactly the same 
topological origin as the baby skyrmion in Skyrme theory, 
but now the topological quantum number is expressed by 
TT2{S^) of the condensate ^, 

q^~^ J e,,d,^^d,^d^x^n. (20) 

This clarifies the topological origin of the non-Abelian 
vortex in two-component BEC. 

The helical vortex will become unstable unless the 
periodicity condition is enforced by hand. But just as in 
Skyrme theory we can make it a stable knot by smoothly 
connecting two periodic ends. In this knot the periodic- 
ity condition is automatically guaranteed, and the very 
twist which causes the instability of the helical vortex 
now ensures the stability of the knot. This is so because 
dynamically the momentum mk along the z-axis created 
by the twist now generates a velocity current and thus 
a net angular momentum which provides the centrifugal 
repulsive force preventing the knot to collapse. 

Furthermore, this dynamical stability of the knot is 
now backed up by the topological stability. This is be- 
cause mathematically the doublet ^, after forming a knot, 
acquires a non-trivial topology 7r3(S'^). And the the knot 
quantum number is given by the Chern-Simon index of 
the velocity potential, 

= J eijkViVjkd^x = mn. (21) 

This is precisely the linking number of two vorticity 
fluxes. As importantly, this is formally identical to the 
knot quantum number in Skyrme theory Hol ITU . This 



assures the topological stability of the knot, because two 
fluxes linked together can not be disconnected by any 
smooth deformation of the fleld configuration. 

We can estimate the energy of the knot, noticing that 
the radius of the lowest energy vortex ring is about four 
times the vortex tube size This suggestes that the 
lightest knot has the energy comparable to the energy 
of the lightest helical vortex in one periodic section with 
k ~ 1/4k. So the lightest knot in ®^Rb is expected to 
have the energy of the order of 3 x 10^'^ eV. 

The existence of a knot in Gross-Pitaevskii theory of 
two-component BEC has been proposed by several au- 
thors 0jlEI3- In this paper we have clarified the physical 
meaning of the knot. Just as the knot in Skyrme theory 
is a twisted magnetic flux ring, this knot is a twisted 
vorticity flux ring. It has a topological quantum number 
given by the Chern-Simon index of the velocity poten- 
tial of the condensate, and enjoys both topological and 
dynamical stability. 

What is remarkable is that this knot is almost iden- 
tical to the knot in the gauge theory of two-component 
BEC that we proposed recently Both are vorticity 
knots whose topology is identical. This implies that we 
have two competing theories of two-component BEC, the 
Gross-Pitaevskii theory and the recently proposed gauge 
theory, which can describe the knot. 

Constructing the knot migh t not be simple, but might 
have already been done 0, uM ■ Identifying it as a vortic- 
ity knot, however, may be a challenging task. A detailed 
discussion on the subject will be published elsewhere p^ . 
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